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X. INTRODUCTION 

For the past ten years, a group of resoarchers-*mathematicians and theoretical engineers, 
centered at, and partially supported by, the Flight Control group at NASA- AMES— have attempted 
to push beyond what was done in the 1960's for linear control theory, and develop effective 
methods for controlling systems whose dynamical equations are fundamentally nonlinear. Our 
applied focus has been the practical problems encountered in designing aircraft and helicop- 
ters, but our methodology— based as it is on fundamental mathematical principles— is adapt- 
able to robotic systems. 

Conversely, we hope that use of the new ideas under development in the computer science 
and AI community will help us use computer technology in a more effective way to handle types 
of control problems — particularly of a "discrete event" nature — that have been difficult to 
include in a differential-equations based methodology. 

Taking a historical view of progress in engineering and engineering-related mathematics , 
the situation becomes clarified. The breakthroughs of the 1960' s in control theory were 
closely linked to the development of computers, which could solve differential equations very 
efficiently. Mathematically, assumptions on linearity worked well because of the nature of 
the engineering problems that needed to be solved, especially in the Apollo Program, where 
the space craft could be treated satisfactorily as point particles, or at worst as rigid 
bodies. In the 1970's we attempted to adapt the mathematical techniques developed in the 
1960 's to the more difficult problems of control of aircraft and helicopters in circumstances 
where the assumptions of linearity of the dynamics could no longer be realistically justified. 
Recently, there has been a change in computer technology — such as LISP logic-based symbolic 
computation and greatly increased potentialities for parallelism — that has not yet been fully 
integrated into the main body of control theory. Further, computer science has achieved 
greater maturity and substance, and I believe that there are great scientific and technological 
possibilities in combining the talents and insights in the two communities. What control 
theory has to offer is a mature, mathematically based overview of a certain class of engineer- 
ing problems, based on concepts of differential equations and dynamics, while the youthful 
vigor of the computer science discipline is generating a lot of energy, but exhibiting the 
need (in my opinion, at least) for more scientific and mathematical direction. 

For the past two years, I have been trying — with George Meyer's advice and support on the 
engineering questions— to push in two directions. First, to understand how the control 
techniques of feedback linearization — developed as a useful control algorithm by Hunt, Meyer 
and Su at NASA-AMES [1,2] — can be integrated into the mainstream of differential geometry and 
extended in the direction of understanding the relation between global and local feedback 
linearization. Second, I have tried to familiarize myself with the LISP and logic-based 
computer technology and algorithms, and help in the job of introducing it into control theory. 
Since the first part of this program is further along — a major mathematical paper is now 
completed [3] and awaits publication — I will describe some of the ideas it contains here, and 
leave my ideas about developing relations between computer science and control theory to 
anotner occasion. 

2. A VIEW OF FEEDBACK CONTROL IN THE CONTEXT OF DIFFERENTIAL EQUATIONS, DIFFERENTIAL GEOMETRY, 
AND LIE THEORY 

A feedback control system can be taken as an underdetermined system of ordinary differen- 
tial equations of the following general form: 

f ( x# if ' u ) * 0 (2.1) 

x E R n ; u € R m 
f is a map: R^ n+m -♦ rP 
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i* * vector of r” describing component* of the system (aircraft, helicopter, spacecraft, 
robot,. ••) that are fixed in value, such as velocities, positions, angular or linear momenta, 
etc. -u- are the control variables, which we must choose in some way to achieve a prescribed 
goal. 


Feedback c control can be described as follows. A feedback map ,or law is a map 

x ♦ F(x) - u (2.2) 

0 R 10 

from an open subset 0 of R n to the control space R m . A trajectory of the feedback 
control law (2.2) is a curve 

t - x ( t) (2.3) 

in R n that satisfies the following ordinary differential equations 

f(x<t>, §f, P(x(t))| - 0 (2.4) 

In engineering practice, we will want to choose the feedback law (2.2) , so that the 
family of trajectories defined by (2.3) and (2.4) will have certain stability, robustness, and 
design properties. (For example, for the latter one might want the trajectory (2.3) to start 
off at time t ■ 0 at a point xq and end up exactly or approximately at a point x A at t-t^.) 
Stabilisation is the property that is best understood mathematically, hence I will use it as 
a touchstone here. 

Much of the work in control theory of the 1960's — which was very successful on both the 
mathematical and practical fronts — was oriented toward linear control systems, i.e., those of 
the forms 


3 ^ - Ax - Bu - 0 (2.5) 

where A and B are constant matrices of appropriate size. Here, it is natural to require 
that the feedback (2.1) preserve this linearity. This can be accomplished by specifying that 
the feedback map ( 2 . 2 ) be of the following form: 

u - Kx ( 2 . 6 ) 

where K is an m*n real matrix. The trajectory equations (2.2) are then of the following 
form: 


|| - (A + BK) (2.7) 

One may then require that these trajectories have a prescribed degree of stability. Because 
(2.7) is a system of differential equations that can be handled with well-Known mathematical 
teenniques, we know that this behavior can be specified by imposing conditions on the eigen- 
values of the n * n matrix 

A+BK (2.8) 

In turn, this "pole-placement" problem can be handled with well-known mathematical techniques 
(matrix Riccati equations or Kronecker pencil theory) that were applied in the 1960's, but that 
of course go back many years in the mathematical literature. 

It is especially interesting that the useful sufficient conditions for stabilization of 
( 2 . 5 ) via linear feedback ( 2 . 6 ), i.e., "pole-placement" in the engineering jargon, involve 
controllability of the control system ( 2 . 5 ) and is a mathematical concept that is — as I showed 
many years ago (4] — essentially differential- geometric and Lie-theoretic in nature. Thus, it 
should be no surprise that the problem of stabilization and feedback control of a more general 
nonlinear system of type (2.1) also involves differential geometry and Lie theory. 

Indeed, the work of Hunt, Meyer, and Su [1,2) (preceeded by work of Krener [5], Brockett 
[ 6 ], Sommer [7), Jakubzyk and Respondek ( 8 )) demonstrate this in a decisive way. Their work 
only dealt with feedback control of a certain class of systems (the feedback linearizable 
ones, with the functions f( , , ) occurring in ( 2 . 1 ) satisfying certain conditions) if the 
trajectory stayed within a small neighborhood R n , whose size could not be specified in 
advance. This posed the question of finding conditions for global feedback equivalence. 

There has been important partial work on this problem by Boothby, Dayawansa, and Elliot [9-11] 
using the tools of differential topology and foliation theory. In my paper [3] I have begun 
to develop ways of applying the Ehresmann-Haef liger [12] theory of pseudogroup cohomology to 
this problem, but there is a long way to go before the results that are useful in practical 
situations will come forth. 
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The mathematical ha art of the mat ho da I have davaXopad in 13] is tha thaory of vactor 
fiaid systems (or distributions) on a manifold and thair equivalence* X will now sketch soma 
of this basic differential-geometric thaory, than raturn to tha control situation* 

3. VICTOR FIELD SYSTEMS AND FEEDBACK EQUIVALENCE 

X will now usa tha formalism "calculus on manifolds," particularly tha thaory of vactor 
fialds (i.e#, first-ordar linaar partial diffarantial oparators) and tha Jacobi-Lie bracket 
[ . ] (i.e., commutator) of such vactor fialds* Saa Xsidori's book (13) for an engineers 
introduction to thasa concepts* 

Lat 8 be a manifold, with V(8) tha space of vactor fields* In terms of coordinates 
(»*) for 8, 1 < i,i < M ■ dim z7 a V€V(Z) is a differential operator of tha following 

forms ” 

V - A 1 (*) $-r (3.1) 

a * 1 

(summation convention in force) 


If 


than 


V* 



[V,V] 



- B 


3A^ \ 

7 ?) 


a 

57 


(3.2) 


(3.3) 


Let F (Z) be the ring of C* , real-valued functions on Z. V(Z) is a module over 
F (Z) , since vector fields can be multiplied by functionss 

(f,V) * fA 1 (3.4) 

a* 

Definition , A vector field system W on Z is a subspace of V(Z) satisfying the follow- 
ing condition : 

fV € W for V € W 
V^V 2 € W for V 1 ,V 2 € W 


i.e., W is a submodule of V(Z). 

Let W be such a vector field system. For z € Z, set 

W(z) * (V(z) : V€ W) (3.5) 

W(z) is a linear subspace of the space of tangent vectors at z. Its dimension is called 
the rank of W at z. W is said to be nonsingular if the rank is constant as z ranges 
over Z. In this paper we will assume that all vector field systems considered have constant 
rank, unless specified otherwise. The concept defined next will play a basic role in this 
work . 


Definition . Let W c V be a vector field system. Set 

C(W): (V € Ws (V,WJ c W) (3.6) 

C(W) is called the Cauchy Characteristic system associated with W. 

Theorem 3.1 * C(W) is another vector field system on Z with the following properties: 

C(W) c W (3.7) 

(C(W) ,C(W) ) c C(W) (3.8) 

i.e,, C(W) is Frobenius integrable as a vector field system 

(C(W) ,W] c w (3.9) 


Proof. Follows from (3.6). 
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Definition , A curve t • *(t) in z is called an orbit ourv # of tha vactor field lyitta If 
if *We following condition ia satisfied! 

Thar a la a vactor flald V 


V 



In W auch that t + s(t) la an orbit curve of V, i.e., if 
|| - V(a(t) ) 


( 3 . 10 ) 


or, in coordlnata terms i 


- A(z (t) ) 


( 3 . 11 ) 


In tnia way, a vactor flald ayatam dafinaa a family of curvaa on Z. It ia thia gaonatric 
proparty that la tha key to tha uaafulnaaa of vactor flald ayatama in control theory. Aa we 
have aaan, control ayatama ara alao defined by familiaa of curvaa, namely aolutiona of tha 
control aquation* i 

f ( x(t) * S» u(t) ) ■ ° (3 * i2) 

x e r" , u € rP 

Sett 

Z m R n * R® 


Z ■ (X,U) 

Wa can than define a vactor field syatem W on Z aa tha amalleat submodule of V(Z) whose 
orbit curvaa are solutions of tha control equation (3.12). 

Let X and X 1 be manifolds. Let W and W' be nonaingular vector field systems on 
X and X', respectively. Let 

a: X - X* 

be a dif feomorphism. 


Definition , a is called an tquivalenoe from the vector field system W to the vector field 
system W 1 if the following condition is satisfied: 

a # (W ( x) ) » W'(a(x)) (3.13) 

for all x € X 

i.e., if a maps an orbit curve of W into an orbit curve of W'. Our problem is to describe 
numbers attached to vector field systems that are invariant under equivalence. We shall cite 
(without proofs) some of the theorems from (3) that do provide such invariants. 

Theorem 3. 2 . Let a: X X v be a dif feomorphism from X to X* that is an equivalence of 
vector field system W to vector field system W*. Let C(W) and C(W*) be the Cauchy 
characteristic systems of W and W*, respectively. Then, “the following condition is 
satisfied: 

a* (C (W) ) - C(W*) (3.14) 

i.e., a is an equivalence between the Cauchy characteristic systems of the given vector field 
systems. 

Definition . For the vector field system W, set: 

W 1 - W + (W,W] (3.15) 

It is called a derived aystem of W. 

Theorem 3.3 . Let a be an isomorphism from W to W*. Then, it is an isomorphism of the 
derived system w* to w' 1 . 
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We can now iterate. Set 


(W 1 ) 1 ■ 


(3.16) 


to define the euooeeeive derived eyeteme of the given vector field syetem W# denoted ae 
*r#W #••• we obtain an increasing filtration of submodules of the module of all vector fields 


on Xt 


W c W 1 c W 2 c ... 


(3.17) 


Theorem 3.4. We havei 


C(W) e C(W X ) c C(W 2 ) c ... 


(3.18) 


In words# this says that the Cauchy characteristics of the derived systems also form an 
ascending# filtered sequence of submodules of the module of all vector fields on X. 

We assume that all the modules (3.17) and (3.18) are of constant rank, set 

r ■ rank W 


c - rank C(W) 
r, - rank w 1 


(3*19) 


rank C(W‘) 


and so on. 


Theorem 3.S . The sequence of integers 


r i r l £ r 2 i * * ' 


( 2 . 20 ) 


c 1 c ! £ c 2 i 

attached to the vector field system W are numerical equivalence invariant*. 

Let us now apply these results to control systems in state space form. 

4. FEEDBACK INVARIANTS FOR CONTROL SYSTEMS IN STATE SPACE FORM 

Let us now specialize the feedback control system to consider those of the following etate 
apaoe form: 

- f(x # u) (4.1) 

y e R n , u € R ra . 

Theorem 4.1 . Let 

Z - R n *R m 

- { ( x , u) : x € R n # u € R m ) 

Let W be the vector field system Z generated by the components of the following vector- 
value3 vector fields on Z : 

s - ! f(x ’ u) k'h\ (4 - 2 > 

Then# the orbit curves on W are precisely the curves t *•> (x(t)#u(t)) that satisfy the 
control equation (4.1) , . 

Theorem 4,2 . Let dx/dt * f(x#u)# and dz/dt » h(z#v) be two feedback control systems with 
tne same number of states and controls. Let 


! f (x,u) 


i_ JL I 

3x 3u | 


and 
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W' - jh(y,v) , I 7 J 

b% the vector field systems assigned to these control systems. Let 

Ti R n+m R n+m (4.3) 


be a C* 


with 


map of the following form: 
T(x,u) «► <y,v) 


y - a(x) 
v - 0 (x , u) 


(4.4) 


(4.5) 


Then T maps the control system {dx/dt * f(x,u)) into the control system 
(dy/dt ■ n(y,v)}, in the sense that it maps solution curves of the first system of ordinary 
differential equations into solution curves of the second, if and only if T is an equiva- 
lence from the vector field system W to the vector field system W'. In particular, the 
integers r , ri , • • . ; c,c^ , . . • assigned by (3.3) to W are invariant~under feedback equivalence. 
In the case of a linear control system, these integers can be computed in terms of the 
aontrollabi l i ty indioee. 

Let us now consider the vector field systems associated wi:h a linear, scalar input, 
control nystem, i.e., one of the following form: 


dx 

It 


Ax + bu 


(4.6) 


x € R n , u € R, b € R n 

Associate with that system the following pairs of vector fields on R n : 

» ■ “k 


(4.7) 


- b i- 
3x 


Let W be the vector field systems on R spanned by these two vector fields and 3/3u. 
Set: 


V A =» Ad (V) ( V Q ) 


Iv ' v i-i> 


(4.8) 


for i > 0 


Theorem 4.3 . The following commutation relations hold among these vector fields on R : 


[V,V i l > V i+i , for i - 0,1,2,... 


(4.9) 


(v i ,v j ] = 0 , 


for i , j * 0,1... 


The derived system is the vector field system generated by 

I tu ' V#V i ; 1 * if * * m 'i ! 


(4.10) 


Theorem 4.4 . If the system (4.6) is controllable, then 




I Tu 


V V 1 



(4.11) 


As I show in [3], Theorem 4.4 
Hunt, Meyer, and Su [1,2] provided 


is the geometric heart of the sufficient conditions that 
in their work on feedback linearization, namely: 
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Theorem 4.5 . Let be vector fields on R n generating a single input, controllable 

control system of the rol lowing form 

- V(x) + UV Q (X) (4.12) 

Suppose the following condition is satisfied: 

The vector field systems 

{ V V l“ v j * (4 - 13) 

are Frobenius integrable for all j . 

Then, the system (4.12) is locally feedback equivalent to a chosen system. 

In the Hunt-Meyer-Su work, the transformation T, which establishes the feedback 
equivalence of (4.12) with a linear system, is obtained as a solution of a system of first 
order, partial differential equations, and we can only prove existence of such linearizing 
transformations locally . A basic question is: 

How to piece together such local feedback equivalences 
to find a global one? 

The answer can be described in terras of cohomology theory 112) . Indeed, this is a typical 
problem of global differential geometry: 

Find conditions for the existence (and computational 
feasibility!) of a global solution of a system of 
partial differential equations when the conditions for 
existence of local solutions are satisfied. 

What complicates the analysis of the conditions for existence of a global solution is 
that the cohomology theory one must use involves an algebraic object — the groupoid of feedback 
automorphisms of the linear control systems — that is infinite dimensional , so that standard 
topological techniques are not very helpful. It is interesting to note that elementary 
particle physicists at the frontiers — in the so-called string theory — are involved with 
mathematical monstrocities that are very similar to these! Work on this question is in 
progress. 
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